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Summary 

A fast coopuMr proip«m. ORID3C. has been 
develop ?d for accurately •enerMing periodic, 
bouadary^conforming. thrce-dimeasional. 
consecutiv^ reflMd comiwtttkMMl grkb applicable 
to axiid litfboaMdilaery gaoaa^^. The 

method U baaed oa ittiag two aaal^ ftmciioas to 
gentfate two-diiiienaioaal ^kis on a number of 
coaxial axisynunetric surfaces poekkmed between the 
cemerbody and the outer radial boundary. These 
boundtf;, fitted grids are of the C*type and are 
chara^erued by quaii'OrthogonaUty and geometric 
periodicity. The buUt*h> nonorthogonal cowtfinate 
stretchings and shearings cause the grid du^ering in 
the regions of interest. The riretdiing parameters are 
pmt of the input to ORID3C. In its |»esait version 
GR1D3C can gemo-ate and store a maximum of four 
consecutively refined threc'dimensional grids. The 
output grid cotMttinates can be cdculated either in 
the Cartesian or in the cylindri^ coordinate system. 


introduction 

When solving nonlinear partial diffwreaitial 
equations governing Huid How. exact boundary 
conditions should be andied on all boundaries. This 
requironent is easily met with a ^id generated by the 
computer program ORID3C, which creates a three- 
dimensional C-type grid conforming to an irregularly 
shaped blade (or wing), hub (or fhselage), ai^ duct 
(or outer radial bound^ or wind tuni^ wdl). 

Three-dimensional boundary-conft^ing grids for 
cascades of Uades (as well as for the single wing- 
body combinatkm) can be gtmerated by a variety of 
methods (ref. I). The method used in the |»esent 
work (ref. 2) has the distinct advantage of being very 
fast and easy to und^and. It is based on two 
conformal mapikng functions (ref. S) and on several 
additional simple uialytical rriatkms ilwt perform 
coordinate stretchings and shrarings. Ilte gmeral 
concept is similar to an rarlicr work on dev^ping 
three-dimensional grids of an O-type (refs, t mid 4). 


*Prm«itly Asristant Profmsor of Amo^Mice Engi- 
neering and Engineering Mechanics. The Univmiity 
of Texas at Austin. 


The CogyutiiiottalFhddMadsMtict Branch trfttse 
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h^fhl in t^aiataig the three-dlmauional p^. 


ApplioibiUly of the Computa* 
Program 

CoBqnttar pragrwi OR1D3C U cigi^Ma of 
genartthig xvy.t or x,$,r coordtottct of a maxhmmi 
of four Gomacutiv^ rdfoad botmdary-confomdng 
grkU for any of tIM foDowii^ UuraiHBaMMioiial 
coaflguratkms: 

(1) ArMtrarUy sh«|>ed wing rakimounted on an 
axisynunetric ftmdage: This wing-body otunbiaa^ 
can be in free air or caitraUy i^tkmed in a wmd 
tunnel having a doubly infinite axisynunetric wall. 

(2) Propeller-type (usually referred to as 
"horizontal axis") windmill rmor having an 
arbitrary number of arlntrarily shaped Madm 
mounted on an Axisynunetric hub 

(3) Helktqtter rotor 

(4) Proprilo' (prop-fan) for airaaft luopulsicm: 
The proprilm can be in (Tee air or omitraUy 
positioned in an axisymmetrk wind tunnel. 

(5) Axial turbomachinery stator 

(6) Axial turbomachlaery rotor 

(7) Ducted fan with rotor diameter sroaUer than 
the duct ^ametm 

(8) Propelim for marine propulsion: The propeller 
can be either free or ducted. 

In the case of a turbomachinmy stator, rotor, or 
ducted fan the outer rmhal boundary represents the 
surface of a duct or shroud. In the remaiidng text this 
surf^ is referred to as a duct. In the case of a free 
propriler, a hrilctqstar rotor in hovm, or a wii^body 
combination the outer radial bmindary represtmts an 
mbitrarily axisymmetricatty shaped surfMC po- 
sitioned off the Made tip. 

It is important to note that both the hub (m- 
centerbody) and the duct (wind tunnri wail <u outm 
radial boundary) can have different axiiymmetric 
shapes. The number of biadm is thcmeticaUy 
arbitrary, althou^ ORID3C has certain limitations, 
which are exf^ined in the fMlowing section. Badi 
blade can be defined by as many as 43 diffment local 
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•irfoU sections. The genernl Made (or vring) can 
incorporate an arbitrvy qwwlw of 

chord lenith, thteknm, twte an^, sweep an^t and 
dihedral an^. ORII^ k eqM^ to 

blades (or wi^) with blunt (or rounded), w e dg e d , 
and cusped trailing and/or ieatUng ed^. The 
trailii^ edge on be t^Mi or doeed and the grid wiU 
conform with the approximate wake shape. 

htote thm in the case of a wing4>ody combination 
or a propeller the loni airfoil shapes at the off-tip 
locations do not have to be spedfM in the input to 
ORID3C. This n>anwke extension of the blade shape 
off the tip will be performed auttunatkatty by the 
code while keqtins the same fig>>tOKAord rmio and 
airfoil shape u qsedfied at the Made (or wii^ tip 
location. 


Limitations of the Computer 
Program GR1D3C 

Computer program GRID3C has several 
limitations. They are the fdlowing: 

(1) All geometric parameters defining the Made 
siMpe must vary smoothly in the spanwiw diredkm. 

(2) The hub and ritroiM are doubly infinite. If the 
actual hub has a finik lotgth, it must be specified as 
having thin “stings’* extending from the end points 
toward the axial infinities. 

(3) The hub and shroud radii must vary mKXMhly 
in the direction of the axis of rotation. 

(4) The blade tip chord length must have a 
nonxero value. 

(S> The maximum number of blades is limited by 
the fact that the minimum local value of the M 4 >*to- 
chord ratio at any spanwise station must be greater 
than approximatdy fh/cVin a0.6S. whne h is the 
arc length spacing between the blades and c is the 
local btauie chord loigth. This limiting ratio is higho* 
if the blade Mctkms we highly staggmd or the 
airfoils have higher camber and rriative thickness. 
The maximum number of Mades NB cm be a pHori 
estimated as 

NBa2v(^) : - r 

\c/ (h/c)mi„ 

where r is radial distance from the axis of rotation, 

(6) The maximum absMute value of the local 

stagger an|^ (angle bwween chcml line and 
meridional plane) fw tyiMcai axial compressor blades 
is limited in GR1D3C to The value of 

idinu can be considwrably larger (K^.«90*) in the 
outer portiems of wind turMne Mwlet or pr^pdlm 
where (h/c) is also very large. In an axial tiu turbine 
blade row. will be lower because of the low 
values of (A/r)n^ and the high camber and relaiivc 
thickness. 

(7) GRID3C <k>w not Mlow for the existence of 
part-span dampers between blades. 


(I) Betides Mades. hub, and shroud, no other 
et^d bOMdaries nay be preeen M the (low (Md. 

A dkadvanti«e of the preeent method U that it U 
not appli cab l e for very tMdt, hi^ily twi^ and 
canbotMl Mades that are v«ry doe^ spaced. 


Qrid-Oenerating Concept 

When analyzing any of the configurations 
mentkMMd in the seakrn on ^H^otbUity, it is 
sufficient to consider a ringlc rotwicmally periodic 
segment of the (low fidd (fig. 1). This s e gme nt is a 
doubly infinite periodic volume stretchi^ fas the 
directioa of the mtis of rotation. The volum hm a 
constant angufau width of 360*/NB. where NB kthe 
total number of Mades. The Mades have artitrary 
spanwise distributioitt of ti^, sw eq > , ^KhedrM, and 
twte ai^es. The bad lArf^ diapes can vary b an 
arMtrwy flmhiai along the bbde mmn. The rotor hub 
and the duct (or shroud) can have different wMtrary 
axisymmetric shapes. Such an arbitrary three- 
dimensiooal phyric^ doraafai k fim disoetized b the 
spanwise di'c^on by a numbw of ccNutial 
axkymmetric nirfhoes that are. in general, spaced 
bet v w cn the hub and the shroud. Each of these 
rarfaces k then dkeretbed. 

The major pr oM em in generating and discretizing 
the axisymmetric surfaces is an accurate 
detemtina^ of the imersection contours between 
the irregular blwte surface and the coaxbl 
axkymmetric surfacm coning the blade. The 
coordinates of the pMnts on these contours are 
defined by fitting cubic qdtates abng the Made ami 
interpobting at the radial stations corresponding to 
each axkymmetric surface (ref. 2). 

Once the shape of the btwseaion contour on a 
partkubr cutting axkymmetric surfm k known, the 
(M^lem becomes one of ttiscretizing a doubly 
connected two-dUnenokmal dmnab (fig- 2).SiKh a 
two-dimenrional doniab k mapped crnifwinally on a 
unit circle (fig. 3) with the anahftic function (ref. 3) 

|~e*(2dsin|S+2c<»/5In(2cos/J)l 


w-iw)~2cosd(l-W) 

The unit circle k “unwr^^ed,’’ with the second 
analytic function (ref. 3) 

k-tanh^y^ 

Resulting is a deformed rhombMdal riiape that is 
th«i mMKMthotonally sheared b the horizontal and 
vertical direakHU to give a rectangular X,Y 
computational dkmwb (fig- 4). A uniform grid in tlw 
X, Y plane that is symmttriatUy spaced with respect 


2 



to the Y axis remaps back into the physical plane as a 
C-type boundary-con forming grid (fig. 5). The actual 
radial coordinates are obtained by Htting cubic 
splines along the C-type grid lines and interpolating 
at a number of axial stations at which the radius oT 
the corresponding axisymmetrk surface is known. 


Computer Program 

General Dcecr^Uea 

The computer program GR1D3C consists of a 
main program (fig. d) and four subroutines. Input to 
GRID3C is accomplished by a separate deck of 
cards, which will be discussed in detail in the 
following section. Parts of the output from ORID3C 
appear in standard computer printout form; the 
major portion of the output (representing 
coordinates of each grid point) is automatically 
written on tapes or disks (units 11, 12. 13, and 14). 

The main program reads all input data and 
determines coordinates of the intersection contours 
created by the axisymmetric computational surfaces 
and the blade surface. This is pCTformed by fitting 
cubic splines (subroutine SPLIF) along the blade 
span and interpolating (subroutine INTPL) at 
desired radial stations (ref. 2). The intersection 
contour on each axisymmetric surface is ex|UTSsed in 
sheared cylindrical coordinates and iteratively 
mapped by using subroutine TRANSF onto a 
deformed unit circle. Subroutine TRANSF 
“unwraps” the circle with the semislit. The resulting 
deform^ rhomboidat shape is then transformed into 
a parallelepipedal shape by using nonortht^onal 
coordinate stretching and shearing in both 
coordinate directions. Subroutine TRANSF maps 
points back into the sheared axisymmetric surface. 
Coordinates of the resulting body-conforming grid 
on each periodic strip are then interpolated on their 
respective axisymmetric surface. Subroutine 
XYZINF suitably shapes grid cells at axial infinity 
cutoff boundaries. XYZINF also generates periodic 
and imaginary grid points (fig. 7). Finally, x,y,z or 
x.9,r coordinates are written on disks or tapes 
numbered unit II, 12, 13, and 14. 

laput 

An example of the entire input to GR1D3C, in card 
form, is shown in figure 8. llie first card is the title 
card on which the user can write as many as 80 
alphanumerical characters (fig. 9) specifying the 
name of the input deck. This text vdll also appear on 
the output listing from GR1D3C. 

The second card contains: 

NB number of blades. The minimum value is 
NB>2. 

NX number of grid cells to be generated on the 
blade surface, that is the number of grid cells 


on the first (ooarte) grid in tl» direction of the 
conmHitatkNwl X axk. ^iggeiMd value is 
NX-20. Users of GRID3C should remem- 
ber that NX must ba an even mmb«r. The 
naximian v^ue fmr NX is 180. 

NY number of grkl cdls to be generated in tlw 
computation Y dknetion, Uuu is, the num- 
ber of C4ayers of grid oeUs around the 
idale on the fbrst (oe«se) grid. Suggested 
rn'iidium value is NY-S. It is advisal^ to 
st>cdfy the vidM of NY tnm the fdlowing 
r..!a‘.ion: 

NY -NX/4 

The maximum value is NY -40. The number 
of grid cells between the tniling edge 
and the downstrmm m^hy (wake cells) is 
automatiaiUy assumed as eqtud to NY. 
Hence, the tmal t,umber of ptd cdls mi any 
C-tayer a 

NXX-NX+2*NY 

NZ number of grid c^s to be gmerated in the 
spanwise direction («>mputatk>nal Z direc 
tion) from the hub to the outer radiid bound- 
ary on the first (coarse) grid. The shape and 
radial distance of the outer boundary surface 
will be qiedfied later in the input. The maxi- 
mum value is NZ-40 

NT nuffibo' of grid cdls to be generated in the 
spanwise direction from hub to tip on the firM 
(coarse) grid. The value of NT sitould be 

NTsNZ-3 

NT is equal to NZ only wh«i r, is equal to the 
radius of the outer boundary, where r, is the 
rotor radius. 

NH number of input points defining hub and sting 
surface. The maximum value for NH is 43. 

ND number of input points defining duct (or outer 
radial boundary) surface. The maximum value 
for ND is 43. ND must be equd to or less than 
NH. 

NP number of input planes. The maximum value is 
NP-43(ng. 10). 

NG number of grids to be goierated. Ii NG- 1, 
only the coarse three-dimensional ^ J consist- 
ing of (NXX*NY*NZ) grid cdU will be gener- 
ated and its x,y,z (or x,9,r) coordinata perma- 
mmtly stored on unit 1 1 . If NG - 2, this coarse 
grid wiU be automatically rcfuied by douNing 
the number of grid cdls in each of three 
computationd directions. Coordinates of this 
refined (second) grid will then be written on 
unit 12. but the coordinates of the coarse 
(first) grid will r«nain stored on unit 11. If 
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NG»3. the Meond grkl will also be refined hi 
the same naiuior at the fint trid. The go* 
ordinates of this fine (third) grid will then be 
writtn on tti^t 13. If NO«4, the sa<ne pro* 
cethtfe will be r^mted, and the coordlnun 
of Uw fourth (^al) grid will be wrh^ oa 
unit 14. Users of OR1D3C should be atwure 
of the computer storage requiienients result* 
ing frmn the SHCcesdve grid refincnmit 
procedure (fig. 11). The present version of 
GRIDK; is capaUe of automatiadly gener- 
ating Goordiiitues fo: as mwiy as four 
(NG«4) refined grids, the fourtii grid hav- 
ing a maximmn of 240x40x 40 ^id cells. 
These values can be changed by duui^ng 
COMMON and DIMENSION statements. 

NC coM^inate system p r e f er en ce tiidcx. If NC» 1 
outpiu coonlhiates will be Cartesian 
If omput coordinates will be cylii^- 
cal (x.M. 

XX X coordinate stretching parmneter (ref. 2). If 
XX* 1.0 ffrid pc^ts on the Mate surface 
will be equidistantly placed with reflect to the 
blMk surface arc length. If XX> 1.0 the grid 
points will be symmetrically clustered 
closm^ to the Imxling edge and the trailing 
edge. The input value for XX must be in the 
region 

I.40>XX2l.0 

YY Y coordinate stretching parameter (ref. 2). If 
YY * 1 .0 the C-type grid layers will be mildly 
clustered towards the blade surface. If 
YY>1.0 the C-type grid layers will be clus- 
tered closely to the Made surface and the 
wake. The suggested value is 

1.05aiYYi0.95 

ZZ Z coordinate stretching parameter (ref. 2). If 
ZZ « 0.0 the axisymmetric computational sur- 
faces will be equidistantly spaced in the axial 
direction. If ^>0.0 these surfaces will be 
symmetrically clustered closer to the hub sur- 
face and the blade tip region . The sug^ted 
value is 

0.75 2ZZ 20.0 

DZ length scaling coefTicient. DZ multiplies all 
output length coordinates generated by 
GRID3C program. 


Input Imsth 
units 

Output Icntth 
units 

DZ 

meters 

meters 

1.000 

feet 

feet 

1.000 

fMt 

meters 

0.3041 

meters 

feet 

3 . 2 m 


RT rotor raitfus, w wing half ^>ao f/. lei^. 

SA sett^ (idtte) ani^ of tlw Mate, teg. $A 
r epre se nt s the constant angk timt will be 
autoaunici^ added to the local stinger (twist) 
a^ tt ea^ qianwte ii^t station when 
genertting a seqiMmce of grids tet a vwiaMe- 
^tdi rotor (or stamt). Each time it is neces- 
sary to diange the valM of SA cmly; the rest of 
the h^Hit should stay the same. 

XE distance betwen the Made trailing edge and 
the down^ream boundary of the flow do- 
main. XE is mipressed in average blade 
diord tenths. 

A number of the ft^iwii^ input ouds (card 4 to 
card min fig. 9) specify hub and duct ^ometry. Each 
these in^t ca^ liMs the following values: 

Xf, X comrdiiuue of the iiqiut pMnt <m the hub or 
AingsurfiKe 

fA r coordinme of tlw iniMit point <ni the hub «r 
sting surface 

x^ X comrthnate of the iiqiut point an Uw outer 
radial boundary 

tfi r eoM’diaate of the inirat point on the outer 
radial boundary , 

Hub and shrtxui surfaces should be defined at least a 
distance of ( 2 «-/NB)r 4 upstream of the rotor. 

Blade geometry is defined on NP imrallel input 
{danes (fig. 10), where the first iniHit plane (N* 1) 
must be entirdy insite the hub a^ the last input 
(dane must be entirdy at beyond the shroud sur- 
face or outer radial boundary sinface. 

The ftdlowing input card (card m + I in fig. 9) con- 
tains gemnetric paranwters specifying the first in|Hit 
plane (N* 1). These parametm are 
ZL X coordinate of the blade leading edge on 
the NOi input plane (fig. 10) 

XL X coordinate of the Mmie leading edge on 
the input plane (fig. 10). Tlw origin 
of the x,9,r coordinate system is 
arbitrarily p(»itioned on the x axis. 
Blade stacking axis corresponds to z 
axis (or r axis). 

YL y coordinate of the blade leading edge on 
the NOi input plane (fig. 10). 

CH chord length (not the projected chord 
length) of the blade on the AA>> input 
plane (fig. 12). Altlwugh the chord 
length is always a positive number, 
the value of CH can be assigned a 
negative sign, thus indimiting if a 
sinite airfoil shape is to be used for the 
blite design. If CH>0 on the 
input Mane, the airfoil shiqw on the 
(N-»-l)di input plaiw is the same as it 
was on the input plane. If CH<0, 
the airfMl shape on the (N-f i)th input 
Mane is different from the airfoil shape 
on the Mh input plane (fig. 13). 
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TW tntle between the blede chord on the 

input pMuie tml the jr axis (f%. 12), dcf. 
TW » positive on Hgitre 12. 

IPNT numhor of umhu points deflnini airfdl 
shape on input plane. The input 
points are given in the clockwise 
direction starting from the trailing- 
edge point. If the airfoil shape is non- 
symmetric. IPNT must be an evoi num- 
(counting the traUing-«ige pdnt 
twice). If the airfoil shape is symiMtric, 
IPNT must be an odd number, and only 
coordinates of points on the loe^ 
surface of the airfoil should be spedfied. 
The value of IPNT cut vary from one 
plane to the next. Maximum number 
is IPNT » 167. 

The following set of input cards (card m 2 to card n 
in Tig. 9) specifies the x’ and input cotxrdinates 
(fig. 12) deflning the airfoil shape tm the NUi input 
plane. The origin of the coordinate system is 
arbitrarily positioned with respeix to the airfoil 
leading e^ (fig. 12). The input coordinates x' and 
y' do not have to be normalized with respect to the 
chord length. Scalings, rotation, and translation of 
the input coordinates will be performed auto- 
matically by the GRID3C program. The coordinate 
axis x' don not have to be parallel with the x 
comdinate axis, but it must be in the TVth input plane. 

All the x' coordinates arc specified in the input 
separately, followed by all the y' cowdinates 
(figs. 8 and 13). The number of x' (or^'') cocmlinates 
that can be specified on each input card is eight (fig. 
9). If the total number of input points (IPNT) is not a 
multiple of eight, the remainder of the last input card 
of a set specifying x' coordinates (as well as the last 
card specifying .v' coordinates) should be filled with 
the values x' - 0 (or »0). 

This concludes all the input data needed for the 
A/th input plane. 

The (/y-f l)ih input plane and the airfoil shape are 
defined on the following set of input cards, starting 
with card n -f I in figure 9. These cards are identically 
formatted as are the cards from m > I to n (fig. 9). If 
CH >0 on card m-f I, the (N-t- 1)<*' input plane will 
be defined by a single input card having the format of 
card m -f 1 . as can be seen in figure 13. 

Output 

Parts of the output of GRID3C will appear in 
printed form and require approximately three pages 
of computer printout (fig. 14). The main result of 
GRID3C are coordinates of three-dimensional body- 
conforming grids (or grid, if NG*I). The 
coordinates of the coarse (first) grid will be 
automatically written on tape (or disk) unit 11. If 
NG = 2 in the input data, coordinates of the second 
(refined) grid will be written on unit 12; and if 
NO ^ the coordinates of the third (fine) grid will be 


written (m unit 13. If NOo4. the coordinates of the 
fourth (final) grid will be written on unit 14. This 
Morage airai^ement provMes ttte tarn wRh an 
oppoHunity to dis|day and analyze Mch grid 
separately. All the WRITE statements arc 
unfomatted and could be found in the midn 
program. 

11w output listing specifies the data contained on 
the firtt m infMit earth as well as each (m-t- l)^eard. 
This serves for easier detecting of mors t^ may 
have been made during the |Mr^ndk» of die iigwt. 

Tlie output from GMD3C continttts with a Hsting 
of the foikwving parameters deffatii^ airfoil shapes 
<m eadi of the K computatkmal axisymmetric 
surfaces (fig. 14). 

XLEAD X coordinate of Made leading edge on Kth 
vurfara 

ZLEAD ;; coordinate of blade leading edge on Kth 
surface 

YLEAD y coordinate of Matte leading edge on Kd> 
surface 

CHORD chord length of Matte (side jMrojectitMi on 
x-x plane) <m Kdi surface 

CHORD - {(yn'-yi^ (JfTE- 

TWIST an^ bctwemi blade chend on K<h surface 
(side iHojection on x-x *xi$ of 

rotation, deg 

TWIST - sin - 1 [(^te ->x4)/CHORDl 

K label (index) number of particular 

computational axisymmetric surface 
generated by ORID3C. The parameter K 
varies from K>2 (hub surface) to 
K*>NZ*NG+2 corresponding to the 
duct or outCT radial boundary surface 
Thus the total number of axisymmetric 
computational surface that will be generated and 
stored on, for example, tape 14 is 

MAXZM-NZ*4+I 

Because subroutine XYZINF generate a layer of 
imaginary points inside the blade and the wake 
(J >MAXYP) e well e one layer (J • 1) of points 
off the periodicity (fig. 7) boundarie (J « 2), the 
total number of C-type layers of points geierated on 
each K<h axisymmetric surface is 

MAXYP-NY«NG+3 

Similarly, GR1D3C geneate grid points on the 
vortex sheet or wake trailing from the blade trailing 
edge (fig. 7). GRID3C also generates two lines (I « 1 
and l*MAXXP) of points off the downstream 
infinity. Hence the total number of points aerated 
on each C-type layer is 
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MAXXP-NXX*NO + 3 

As a conscquoice, tape 1 1 will omiiun coordinates 
of (NXX 3)*<NY -I- 3>*(NZ t) frid points, tape 12 
will contain coordinates of (NXX*2-t-3)* 
(NY*2-f 3HNZ*24-I) grid pomts, tape 13 will 
contain co^inat« of (NXX«3*f 3)*(NY*3 + 3)« 
(NZO+I) grid points, and tape 14 wiU contain 
coordinates of (NXX*4+3HNY*4+3)*(NZ*4+I) 
grid points. These coordinates are written on corre- 
sponding tapes in a sii^le-precistcm unformatted 
fturm. 

To assure the uwr of GRID3C that eMh grid was 
stored mi its corresponding tape or disk, the 
following output control parametm are printed after 
grid was generated on ea^ K<h surfa^: 

NT APE tape (or disk) unit number on which the 
appropriate v^d from K<i> surf^e was 
stored 

MAXXP. total number of grid points (MAXXP* 
MAXYP MAXYP) extracted frwn the finest grid 
on Kt>* surface and stored on tape «^th 
the unit numbm* NT APE 

El«nents of a three-dimensional C-type poiodk 
grid generated by GRID3C code for a simile wing- 
body combination (fig. 8) are shown in figure 13. 

Nme that the physical wii^ is rectangular. Its 
extension in the off-tip rqpon is automatically 
generated by GRID3C while keeping the gap-to- 
chord ratio the sane as at the wii^ tip location. For 
clarity reasons only one of the axisymmetric (in this 


case cyUndrkal) surfMes is shown, and the inter- 
sacfitm cMtoun other surfaom with the wing 
surface are indicated. Figures 16 and 17 represent the 
dtecretoed blade surface and the hub surface of an 
^ght blmle advanrad NASA prop-fan configuration. 
N(He the fact that in this case the hub has a variable 
radius. 
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Figure 2. - Cascade of airfoils (with wakes) on i^,R,6) computational axisymmetric surfaces. Cascade of 
semi-infinite slits and periodic boundaries are indicated. 
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Figure 3 - Intermediate or "circle” plane obtained by conformally mapping a cascade of slits ( ) 

and a cascade of airfoils ( ) from figure 2. 
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F-itturc 5. C'lcncral appearance of a remapped C-type grid on an axisymmctric surface. 
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Figure 6. -OloW How chart for GR1D3C coniiiuter code. 
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Fifure 7. -An example of a complete C-type two-dimensional grid on each of the axisymmetric 
computational surfaces. Note the periodic layer (J > 1 ) of points and the imaginary iayw (J ■ MAXYP) 
of points. Also notice the columns (1-1 and 1-MAXXP) of imaginary poims at the downstream 
boundary. In actual flow cakulatkms, grid point 1 - IMID is omitted. Hence, its location is arbitrary. 
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KECTMieULAIt NACAttll UXN6 mDHWtalTO ON CYL. FttSELAOE : ll^T EM CONIOSCI 
N|: 2NX= 2«MY- »NZ: BNT: SNHsltNO>tailP:2tN«:ENC:l 

XXs I.IBIYY* l.MtZZ: I.SMBZc l.MMT* S.fMSAx O.ltOXEs l.Stt 

-Y.YfatM 0.1590M -Y.YMiM A. •••§•• 

-7.50etA0 •.IStatt -7.5AAMt A.MMM 

•x.tttaai a.isaaaa -s.aaaaaa a.aeaaaa 

-a.aaaaaa a.isaaaa -a.Maaaa a.aaaata 

a.aaaaaa a.isaaaa a.aaaaaa a.aaaoaa 

a.zaaaaa a.isaaaa a.zaaaaa a.aaaaaa 

a.seaaea a.isaaaa a.saaaaa a.aaaaaa 

i.aaaaaa a.isaaaa i.aaaaaa a.aaaaaa 

s.aaaaaa a.isaaaa s.aaaaaa a.aaaaaa 

12 -saaaaa a.isaaaa iz.saaaaa a.aaaaaa 

a.asaaa a.aaaaa a.aaaaa -a.saaaa a.aaaaa a.aaaaa a.aaaaa i.aaaaa a.aaaaa aa 
i.aaaaaaa a.aaiiasa a.aaaaaaa a.aaama a.aaszaaa a.aTsaiaa a.aiaissa a.aiiaraa 

a.aaaaaza a.aaazsaa a.assaoza o.aaaoMa a.assasaa a.azaTiaa a.aizsaia a.aaaaaza 

a.ssaaasa a.saaaaza a.asiaaaa a.os2a29a a.aizaasc o.7aiss<>o i.xaaaaaa a.7a5259a 

a.72oaasa a.aaaassa a.aa7a2t9 a.aaaaaaa a.aizzaaa e.sasst.^ ’ a.ssaasza a.S2sa97a 

a.aassiaa a.aaaasia a.asaaa7s a.aa7ssaa t.s7aaa2a a.ssaxaaa a.szsziaa a.zaasisa 

a.2707asa a.zasaaaa a.222isaa a.iaasasa a.iiazaaa a.isaziaa a.isasaaa a.izziara 

a.iaaiaia a.aaisaaa a.aiazasa a.aaassta a.ass7ua a.aaasiaa a.asaaa7a a.asaaaaa 

o.92a7a7a a.aiaa7aa a.aiaaasa a.aiii7sa a.aa7aaza a.aasasaa a.aazaaaa a.aaaaaaa 

a.aaaaaat-a.aaizaaa-a.aaisaza-a.aaiassa-a.aazssaa-a.aazaaaa-a.aassssa-a.aasfaaa 
-a.aAaaasa-a.aassaaa-a.aaaaaaa-a.aa75aaa-a.aaa8asa-a.aiaiaia-a.aii72aa-a.aisaa2a 
-a.ois22sa-A.ai7i78a-a.ait2aia-a.a2iasaa-a.a257aia-a.a2aisaa~a.a2aaiia-a.a3ii2sa 
-o.aiH7so-a.aia23aa'a.a387aaa-a.aai3aia-a.aas7ssa-a.aaaisia-a.aaasa2a-a.asa5isa 
-a.es2a72a-8.a5a2saa-a.ass7aia-a.a57aaaa-a.as8i23a-a.asas79a-a.asa337a-a.asaaaaa 
-e.esai20D-a.es8838a-a.asaaa3a-a.osaaaaa-a.a95S59a-a.as39a3a-a.as2aMa-a.aaa87aa 
- a. oaisaoo-e.aasaaao-a. 8929888-8.839781 a-a. 8378198-8. a3a22sa-8.83ia2S8-e.82Baas8 
-8. 0239810-8. 8232898-8.8285378-8.8178798-8.8151358-8.8122828-8.0883918 0.0888888 
0.18088 8.88888 8.88888 8.S8888 0.08888 8.00800 8.08800 1.88000 0.08808 89 

8.23808 8.88888 8.80088 0.30080 8.88088 0.00088 8.08000 1.08080 8.88000 89 

8.38080 8.80880 8.88888 0.58809 8.88888 0.88000 0.00009 1.00008 0.08000 89 

0.73880 8.88888 8.08808 0.58080 0.08888 8.88088 8.88000 1.00008 0.88008 89 

1.88888 8.88888 8.88888 8.S8888 8.88888 8.88888 8.88888 1.88888 8.88088 89 

1.58088 8.88908 0.88080 O.SOOOO 0.00808 8.08888 8.88808 1.88888 8.00000 (9 

1.75890 8.00008 0.00008 0.50000 8.88800 8.08800 8.80000 1.08888 0.80080 99 

2.00808 0.00808 8.00888 8.50808 8.80888 8.00888 0.88888 1.00080 8.80800 99 

2.23808 0.00800 0.00000 O.SOOOO 8.80800 8.88800 0.08880 1.88808 0.80088 99 

2.50080 0.88880 8.00808 8.58000 8.08800 0.88888 8.80889 1.88800 0.08008 99 

2.75000 0.00008 0.80000 0.50008 8.00888 0.08808 0.80000 1.08000 8.00000 99 

3.80008 0.00000 0.00000 8.50880 8.00000 8.00000 0.89008 1.80008 8.08000 99 

3.18000 8.00000 8.00000 O.SOOOO 0.80000 0.08000 0.88800 1.00808 0.88800 99 

3.25080 0.90000 0.80000 8.50800 8.80000 8.88088 0.08808 1.80808 8.88088 99 

3 50000 0.00000 0.00080 8.58080 8.88000 8.88808 8.00808 1.88888 0.80008 99 

9.08808 0.00008 0.08890 0.50088 8.80888 8.88888 8.00008 1.80088 0.80008 99 

9.50080 0.00008 8.00000 8.50800 8.00008 8.80800 8.80800 1.80808 0.00008 99 

5.50000 0.08000 0.00008 0.50880 8.08080 0.88888 8.00088 1.08080 0.80088 99 

9.30000 0.00000 0.80800 8.58008 8.88000 0.88088 8.08008 1.80808 8.00008 99 

8.50000 0.00000 0.00000 0.58088 8.00080 0.88888 8.80880 1.88808 8.80808 99 


Figure 8. -Complete input for a geometry consisting of a douUy infinite circular cylinder and a 
midmounted, rectangular, unswept wing composed of NACA 0012 airfoil sections. The wing is at zero 
angle of attack. Note that only the lower wing surface points are specified because the airfoil is 
symmetric. 
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Figure 9. - General input format for GRID3C computer code. 
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Figure 1 1. - Computer requirements and GRID3C performance charts as evaluated on an IBM 370/3033 
computer. 
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Figure 12. -Relation between the global coordinate system and the input Qc',y',z) coordinate 

system. 
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H|s MNX« 24HY> 4N2* 4NTi 4NN«l«HD*llNP«IIN«sSMC«l 
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I.mORI I.SII1M 
l.MMM l.OOMM 
I.IMOM S.44MM 
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1.1141411 
1.2134141 
1.31*4711 
0.4144211 
0.413*22* 

0. 1124170 

1. *44*44* 

0.7447*10 
0.11*174* 

0.174313* 
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9.01*4*40 


(MID3C3 

9. 1*9X1* 1.1*1 
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1 . ****** 
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Figure 13. An imput example describing geometry of a forty-blade ducted rotor with a linear spanwise 
blade chord variation. The blade is composed of two airfoil shapes. 
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Figure 14. - Output from GRID3C corresponding to the input given in figure 13. 
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Figure IS. - Key elements of a three-umensional C-type grid generated by GRID3C code with the input 
given in figure 8. Note the grid deterioration in the ui^tream region. This can be improved by mo^fying 
parts of the main program and the subroutine XYZINF. 


22 


f 


0RU3NAL PA'JaE 
OF POOR QUALITY. 



Figure 16. - Paris of a three-dimensional C-type grid generated by GR1D3C code for an advanced eight- 
blade NASA prop-fan. Note significant spanwise variations of blade sweep angle and chord length. 
Note also that the hub surface is axisymmetric. 
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